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Abstract 

The second reference state of the open XXZ spin chain with non-diagonal boundary 
terms is studied. The associated Bethe states exactly yield the second set of eigenvalues 
proposed recently by functional Bethe Ansatz. In the quasi-classical limit, two sets of 
Bethe states give the complete eigenstates of the associated Gaudin model. 
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1 Introduction 



The open XXZ quantum spin chain has played a fundamental role in the study of quan- 
tum integrable systems with various boundary interactions, which appeared in statistical 
mechanics, condensed matter and quantum field theory. Although the special case of diag- 
onal boundary terms was solved long ago [H El E] , Bethe Ansatz solutions for non-diagonal 
boundary terms where the boundary parameters obey some constraints have been proposed 
only recently by various approaches [U El El CO El E] ■ It is found that in order to obtain the 
complete spectrum of the model two sets of Bethe Ansatz equations and consequently two 
sets of eigenvalues are needed [101 E], in contrast with the diagonal boundary case [3]. This 
suggests that in the framework of algebraic Bethe Ansatz there should exist two reference 
states (or pseudo-vacuum states) corresponding to the two sets of Bethe Ansatz equations 
and eigenvalues. However, to our knowledge only one reference has been constructed so far 
[HI E]- Moreover, the explicit expressions of the complete reference states and associated 
Bethe states are of great importance for investigating correlation functions of the model |llj. 

In this letter we study the second reference state of the open XXZ spin chain with 
non-diagonal boundary terms and construct the complete eigenstates of the model in the 
framework of algebraic Bethe Ansatz. In the quasi-classical limit, they give the complete 
eigenstates of the associated Gaudin model. 

2 The inhomogeneous spin-^ XXZ open chain 

Throughout, V denotes a two-dimensional linear space and o~ ± , a z are the usual Pauli ma- 
trices which realize the spin-| representation of the Lie algebra sl(2) on V. The spin-- XXZ 
chain can be constructed from the well-known six- vertex model R-matrix R(u) G End(V (g> V) 
[IT] given by 



R{u) 



( 1 \ 

b(u) c{u) 
c{u) b{u) 

\ 1 J 



(2.1) 



The coefficient functions read: b(u) = . s , in " N , c(u) = . s , in v , . Here we assume ri be- 
ing a generic complex number. The R-matrix satisfies the quantum Yang-Baxter equation 
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(QYBE), 



R 12 {u x - u 2 )R xz {ui - u 3 )R 23 (u 2 - u 3 ) = R 23 {u 2 - «3)-Ri 3 («i - u 3 )R 12 (ui - u 2 ), (2.2) 

and the unitarity, crossing- unitarity and quasi-classical properties [6] . We adopt the standard 
notations: for any matrix A G End(V) , Aj is an embedding operator in the tensor space 
V ® V ® • • -, which acts as A on the j-th space and as identity on the other factor spaces; 
Rij(u) is an embedding operator of R- matrix in the tensor space, which acts as identity on 
the factor spaces except for the z-th and j-th ones. 

One introduces the "row-to-row" monodromy matrix T(u), which is a 2 x 2 matrix with 
elements being operators acting on V® N , where N = 2M (M being a positive integer), 

T (u) = Rqi(u + z^Rq^u + z 2 ) ■ ■ • R 0N (u + z N ). (2.3) 

Here {zj\j = 1, ■ • • , A^} are arbitrary free complex parameters which are usually called inho- 
mogeneous parameters. 

Integrable open chain can be constructed as follows [3]. Let us introduce a pair of K- 
matrices K~(u) and K + (u). The former satisfies the reflection equation (RE) 

#12(141 - u 2 )Ki(ui)R 21 (ui + u 2 )K 2 (u 2 ) 

= Kziu^Ruim + u 2 )K{( Ul )R 21 ( Ul - us), (2.4) 

and the latter satisfies the dual RE 

R 12 (u 2 - ui)Kf(u 1 )R2i{-ui -u 2 - 2r))K£(u 2 ) 

= Kt{u 2 )R 12 {- Ul -u 2 - 2r 1 )K+(u 1 )R 21 (u 2 - u x ). (2.5) 

For open spin-chains, instead of the standard "row-to-row" monodromy matrix T(u) (\2.3h . 
one needs to introduce the "double-row" monodromy matrix T(it) = T(u)K~(u)T~ l (—u). 
Then the double-row transfer matrix is given by 

t(u) =tr(K + (u)T(u)). (2.6) 

The QYBE and (dual) REs lead to that the transfer matrices with different spectral param- 
eters commute with each other [3]: [t(u),t(v)] = 0. This ensures the integrability of the 
inhomogeneous spin-^ XXZ chain with open boundary. 
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In this paper, we will consider a generic K-matrix K (u) which is a generic solution to 
the RE (j2.4p associated the six-vertex model R-matrix [121 D2] 



^M=(S) (27) 

The coefficient functions are 

2 cos(Ai - A 2 ) - cos(Ai + A 2 + 2£)e~ 2iu 



4 sin(Ai + £ + m) sin(A 2 + £ + u) 
-ism(2u)e' i{Xl+X2) e- iu 
2 sin(Ai + f + u) sin(A 2 + f + u) ' 



^ 2 _ 2sin(2w)e i(Al+A2) e- iu 

x ^ ~ 2sin(Ai + £ + w) sin(A 2 + £ + «)' 

2 _ 2 cos(Ai - A 2 )e- 2 ^ - cos(Ai + A 2 + 2Q 

fc2W " 4sin(A 1 + e + n)sin(A 2 + ^ + n) " 

At the same time, we introduce the corresponding dual K-matrix K + (u) which is a generic 
solution to the dual reflection equation (12. 5p with a particular choice of the free boundary 
parameters with respect to K~(u): 

k+ (u) = ( ) . (2.9) 



The matrix elements are 



k+ i,. _ 2 cos(Ai - A 2 )e-^ - cos(Ai + A 2 + 2Qe 2m+ ^ 

4 sin(Ai + £ — u — rf) sin(A 2 + £ — u — rf) 
, ,i, , isin(2M + 2r7)e^ (Al+A2) e iu ~ f,? 

fc" l " 2 (M) = ; - = ; 1 — = 



2 sin(Ai + £ — u — 77) sin(A 2 + £ — w — 77) 



k+2 , ^ _ -i sin(2w + 2r/)e i(Al+A2) e iu+ir ? 

1 2 sin(Ai + £ — -u — 77) sin(A 2 + £ — u — rj) ' 

x 2 costAx - A 2 )e 2 ^ - cos^ + A 2 + 2§e"^ 

K 2 (w) = : = : : = : — . (2. 10) 

4 sin(Ai + £ — u — rf) sin(A 2 + £ — u — rj) 

The K-matrices depend on four free boundary parameters {Ai, A 2 , £, £} which obey the 
constrain conditions in [H [51 E]- It is very convenient to introduce a vector A = Y^t=i ^k e k 
associated with the boundary parameters {Aj}, where {q, i — 1,2} form the orthonormal 
basis of V such that (e*, €j) = <5y. 
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3 Vertex-face correspondence 

Let us briefly review the face-type R-matrix associated with the six-vertex model. Set 



i 



ti — e, e = \ J2l=i e k, for i — 1) 2. For a generic m G V, define 

mi = (m,€i), rriij = mi - rrij = (m, €i - ej) , i,j = 1,2. (3.1) 

Let R(u, m) G End(V^ ® V) be the R-matrix of the six-vertex SOS model, which is trigono- 
metric limit of the eight-vertex SOS model jH] given by 

2 2 

R(u, m) = J2 R-u(% m)E ll ®E ii +J2 m)E ii ®E ]j +R{;{u, m)E Jt ®E ij ) , (3.2) 

where i£y is the matrix with elements (Eij) l k = SjkSu. The coefficient functions are 

,-,•/,. sin w sin(m,-i — r?) 
J$(«,A =1, flg- «,A = \ ? /; i^j, 3.3 

J sin(« + n) sin(mjj-) 

•i sinr/sin^ + m^) 

Rij{u,m) = — , 1 ± 3, (3.4) 

J sm(it + 77) sin (m^ J 

and rriij is defined in (13.11) . The R-matrix satisfies the dynamical (modified) quantum Yang- 
Baxter equation (or star-triangle equation) [T4] . 

Define the following functions: 9^(u) = e~ lu , 6^ 2 \u) = 1. Let us introduce two inter- 
twiners which are 2-component column vectors 4> m m _ V j(u) labelled by 1, 2. The /c-th element 

of <f>m,m-nj( u ) is S iven b y 



0W ,( u ) = 0(*)( u + 2m,-). (3.5) 

Explicitly, 

/ e -t(u+2mi) \ / e -i(«+2m 2 ) \ 

m ,m-r,i(«) = ( l J > TO ,m-i^( W ) = ( j ) ■ ( 3 ' 6 ) 

Obviously, the two intertwiner vectors <t > m m-m{ u ) are linearly independent for a generic 
m G V. 

Using the intertwiner vectors, one can derive the following face-vertex correspondence 
relation [5J 

#12(^1 - u 2 )4> m,m—rii (Til) ® 

= Yl R ( Ul - M2 ' m )^m- v l,m- v (i+k)( U l) ® 0m,m-r,f( U 2)- (3.7) 
fe,/ 



Then the QYBE ( 12.21) of the vertex-type R-matrix R(u) is equivalent to the dynamical Yang- 
Baxter equation of the SOS R-matrix R(u, m). For a generic m, we can introduce other types 
of intertwiners 0, satisfying the conditions, 

2 2 

SJn-wifa) 0m,m-„*( u ) = <W XI ^m+v[i,m( U ) fi,i),m( U ) = V ( 3 ' 8 ) 
fc=l fc=l 

One may verify that the K-matrices K ± (u) given by (12. 7p and (12. 9p can be expressed in 
terms of the intertwiners and diagonal matrices )C(X\u) and JC(X\u) as follows 

K-(u)t = £>JVj>, A- v -W^(Ak)^?A^(-«), (3.9) 
2f+(ti) t « = E^A-nfM^^I^^CM), *-*(«)■ ( 3 - 10 ) 



Here the two diagonal matrices JC(X\u) and JC(X\u) are given by 

w,i N_ n . / r / \ I \ , a , n n . ,sin(Ai + f-u) sin(A 2 + £-w) , . 
/C A u) = Diag(fc(A u)i, fc A u 2 ) = Diag(— . ), (3.11 

sm(Ai + 4 + it) sm(A 2 + £ + it) 

£(A|it) = Diag(jfe(A|u)i, k(X\u) 2 ) 

sin(Ai2-?y) sin(Ai+^ +u+rj) sin(Ai 2 +?7) sin(X 2 +£+u+r]) , /0 1oN 

= Diagt— — — t t r-, — — : r-77 t H- (3.12) 

sin Ai2 sm(Ai + £ — it — rj) sm A i2 sm(A 2 + £ — it — 17) 

Although the K-matrices if^it) given by (12. 7p and (12.91) are generally non-diagonal (in the 
vertex form), after the face- vertex transformations (13.91) and (I3.10p . the face type counter- 
parts K.(X\u) and JC(X\u) simultaneously become diagonal. This fact enables us to apply the 
generalized algebraic Bethe ansatz method developed in [T5] for SOS type integrable models 
to diagonalize the transfer matrices r{u) (12. 6p . 

The decomposition of K + (u) (I3.10p and the diagonal property (13.121) lead to the recasting 
of the transfer matrix r(u) (12. 6p in the following face type form 

r(«) = tr(K + {u)Y(u)) = £ K(X\uXT(X\u); = £ ~k{X\u) ,T {X\u)». (3.13) 

Here we have introduced the face- type double- row monodromy matrix T(m\u) 

7{m\uX = ^(a-^a W nu)l^ m . v ,{-u). (3.14) 

hi 

This face-type double-row monodromy matrix can be expressed in terms of the face type 
R-matrix R(m\u) (13. 2p and K-matrix K.(X\u) (13. lip (for the details, see equation (4.19) of 
[15]). 
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As in [6], let us introduce operators: 

A(m\u) =A {1) (m\u) =T(m\u)\, B(m\u) = T H^) 2 C r m \ u \ = T (™\ u )i (3 15 ) 

sm(m 12 ) sm(m 2 i) 

V(m\u) =V^(m\u) = sin ( mi2 + r » {T(m\u) 2 2 -R(2u, m + r)l)l\A {1) (m\u)\ . (3.16) 

sin(mi 2 ) "'" 

We remark that the transfer matrix r(u) (I2.6P can be expressed in terms of the operators A^ 1 ' 
and £>W. It was found in [B] that one can construct a reference state, denoted by |fiW(A)), 

|fiW(A)) = A _ (7V -l)^i,A-iV^(-^)®^A-(JV-2) ?? l,A-(7V-lhi(~ 2; 2)---®0A,A-r ) i(-^) 5 ( 3 - 17 ) 

in the sense that the state is common eigenstate of the operators A^ and T>^ and is 
annihilated by C (c.f. (14.61) ). The associated Bethe states can be constructed by applying the 
"creation operator" B on the corresponding reference state 

K, • • • , = B(X - 2 V l\v 1 )B(X - 47711^) • • -B(X - 2Mr ] i\v M )\rt 1 \\)). (3.18) 

If the parameters {vk} satisfy the following Bethe Ansatz equations, 

sin(A 2 + £ + v a ) sin(A 2 + £ - v a ) sin(Ai + £ + v a ) sin(Ai + £ - v a ) 
sm(X 2 +^+v a +rj) sm(X 2 +^-v a -rj) sin(Ai+^+v Q +^) sin(Xi+£-v a -r]) 

_ tt sin(t; a + v k + 2rj) sin(v a - v k + y) 
1£ sin(v a + v k ) sm(v a - v k - rj) 

nsm(v a + z k ) sm{v a - z k ) 
— — — - T , a = l,---,M, 3.19 
fc=1 sra(v a + z k + rj) sm(v a - z k + rj) 

the Bethe state \v±, ■ • ■ , vm)^ becomes the eigenstate of the transfer matrix with eigenvalue 
given by [6] 

A w (u) 



+ 



sin(A 2 + £ — u) sin(Ai + £ + u) sin(Ai + £ — u) sin(2w + 2rj) 
sin(A 2 + £ — u — rj) sin(Ai + £ — u — rj) sin(Ai + £ + u) sin(2-u + rj) 

rr sin(w + Vk) sin(w — Vk — rj) 
J-_-L s in(tt + i> k + 77) sin(w — v k ) 

sin(A 2 + i + u + rj) sin(Ai + £ + u + rj) sin(A 2 + £ — w — 77) sin 2u 
sin(A 2 + £ — u — rj) sin(Ai + £ + u) sin(A 2 + £ + u) sin(2-u + rj) 

n' f sin(w + v k + 2rj) sin(u — v k + rj) 



sin(« + v k + f]) sin(it — i> fc ) 

TT ™( u + -*><) (3.20) 

sin(w + Zk + rj) sin(w — z k + rj) 



k=l 
2M 

X 

fc 
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4 Second reference state and associated Bethe states 



Let us introduce the second reference state |f^ 2 )(A)), 

|fi( 2 )(A)) = A _ ( tf_i)^ )A _tf^(-*l) ® A _ (jV _ 2 ^ (4.1) 

and the associated operators A^ and T>^ which are linear combinations of {T(m\u)l}, 

A^{m\u) = ^^±A {T{m\u)\ - R(2u, m + r}2)HV®{m\u)}, (4.2) 
V (2 \m\u) = T( y m\u) 2 2 . (4.3) 

Using the technique developed in [15], after tedious calculations, we find that the state 
|Q( 2 )(A)) given by (14. ip is exactly the reference state in the following sense, 

sin 2u sin(A2 + £ + u + rf) sin(Ai + £ — u — r/) 



A^\X-Nr]2\u)\n (2 \X)) 



sin(2ii + rf) sin(A2 + £ + u) sin(Ai + £ + u) 



N 



X 



TT + |nC0 (A))j (4.4) 

^ sin(-u + z k + 77) sm(w — ^ + 77) J 

P( 2 )(A - A^^O^A)) = SH ^li~ U \ \n^(X)}, (4.5) 

sm(A 2 + 4 + u) 

B{\- Nr]2\u)\Q {2) {\)) = 0, (4.6) 

C(A- iV772|u)|fi( 2 )(A)) ^ 0. (4.7) 

Then the second set of Bethe states can be constructed by applying the "creation operator" 
C on the reference state |fi (2) (A)> (c.f. (13TTS]) ) 

• • • , v M ) {2) = C(X - 2r ] 2\v 1 )C(\ - Ar]2\v 2 ) • • ■ C(A - 2Mr/2|t; M )|fi( 2 )(A)). (4.8) 

One may check that the transfer matrix r(u) (12. 6p is a linear combination of the operators 
A^ and 

sin(Ai + £ + u + rf) (2) 
sm(Ai + £ — U — r)) 

| sin(A 1 + ^-n)sin(A 2 + ^+n)sin(2n + 2r/) r (2) (A^) (4 9) 
sin(Ai + £ — u — rf) sin(A2 + £ — m — rf) sin(2?j + rf) 

Carrying out the generalized Bethe Ansatz [51 EE], we finally find that if the parameters {v^} 
satisfy the second Bethe Ansatz equations (comparing with the first ones (13.191) ). 

sin(Ai + £ + v a ) sin(Ai + £ - v a ) sin(A 2 + £ + v a ) sin(A 2 + £ — v a ) 
sin(Ai+^+'u Q ,+77) sin(Ai+^-77 Q ,-77) sm(X 2 +£+v a +r]) sm(X 2 +£—v a —r]) 
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_ -ft sm(v a + v k + 2rj) sin(v a - v k + rj) 
±± sin(v a + v k ) sin(v a - v k - rj) 

M h^ M ,„-„) a = 1 ,..., M , (4 . 10 ) 

£J sm(v a + Zk + rj) sm(v a -z k + rj) 

the Bethe states \v±, ■ ■ ■ ,v M )^ yield the second set of the eigenstates of the transfer matrix 
with the eigenvalues, 

sin(2-u + 2rj) sin(Ai + £ — u) sin(A 2 + £ + u) sin(A 2 + £ — w) 



A (2), 



U 



sin(2w + rj) sin(Ai + £ — u — 77) sin(A2 + £ — u — rj) sin(A2 + £ + u) 



rr sin(w + Vk) sin(u — Vk — rj) 
sin(u + Vk + rj) sin(u — Vfj) 

sin(2w) sin(Ai + £ + u + rj) sin(A2 + £ + u + 77) sin(Ai + £ — it — rj) 
sin(2-u + 77) sin(Ai + £ — 77 — 77) sin(A 2 + £ + u) sin(Ai + £ + u) 



n f sin(w + v k + 2?]) sin(u — Vk + 77) 
fc=i 

2M 



X 

fc 



sin(w + v k + 77) sin(?i — 7j fc ) 
sin(w + Zk) sin(w — 2%) 



sin(w + £ fe + 77) sin(w - z k + rj) 



Note that the normalizations adopted in this paper for the R- and K-matrices are different 
from those used in [SJ. After rescaling an overall factor sin(Ai + £ — u — rj) sin(A 2 + £ — u — 
77) sin(Ai+£+u) sin(A 2 +£+7i) Y\k=\ ^.(u+Zk+rj) sm(u—Zk+r]) and setting all inhomogeneous 
parameters = 0, our two eigenvalues {A^(u)} recover those in [9j. Therefore two sets 
Bethe states {\vi, ■ ■ ■ ,vm)^} (13.181) and (14. 8 j) together constitute the complete eigenstates 
of the transfer matrix r(u) (12. 6p . 

5 Results for the Gaudin model 

In order to study the associated Gaudin model, we need further restrict £ = £ [BJ. Follow- 
ing [TBI [T7] one can introduce the corresponding Gaudin operators {Hj} by expanding the 
double-row transfer matrix r(u) (12. 6p at the point it = Zj around 77 = 0: 

d 

r(zj) = id + r ] H j + 0(r ] 2 ), with Hj = — r(^)| r/=0 , j = l,---,N, (5.1) 



where 



2M 1 f (x 2 a z - 1 

^ = r,(%) + E sin( ,. _ Zk) \°t°J + + «»(*,■ - *) JL i— 



where tyu) = ^{i? i (^)}| J?=0 ^(^), j = l,---,N, with = tr {Jr o + (u)i2oi(2u)i2o;(0)}. 

The commutativity of the transfer matrices {t(^)} for a generic 77 implies [Hj,Hk] = 0, 
for z, j = 1, • • • , N. Thus the Gaudin system defined by (15.21) is integrable. Moreover the 
relation (15 .ip between {Hj} and {r(zj)} enable us to extract the eigenstates of the Gaudin 
operators and the corresponding eigenvalues from the results obtained in last section. 
Let us introduce states |flW(A)), 

|OW(A)) = ( e 1 U---®(^ e 1 , 2 = 1,2. (5.3) 

These states can be obtained from the reference states |0^(A)) by taking the limit: |0^(A)) = 
lim^o |^ (A)). Let us introduce a matrix C{u) G End(V^) associated with the intertwiner 
vector (f) (13. 5p 

e -i(u+2X 1 ) e -i{u+2\ 2 ) 



C(u)=^ 1 ^ j, (5.4) 

and the corresponding gauged Pauli operator a ± {u) 6 End(F): a ± {u) = C{u)a ± C{u)~ 1 . 
Then we can construct states ^'(xi, ■ • • , %): 

M / 2M 



,t,(D / x tt / f sm(Ai + i - x a ) sm i a - z k + Ai 

Z-i V I sm(Ai + 4 + x a ) sm(x a - z k ) 

- k (-z k ))p 1 \X)), (5.5) 



a=l \ k =l 

sin(A 2 + £ - x Q ) sin(x Q , + ^ - A i2 ) 
sin(A 2 + £ + x a ) sm(x a + z k ) 

M / 2M 



' M ~i V £^ I sin(A 2 + $ + x a ) sin(x Q - 

fc +(-^) |^ (2) (A)>. (5.6) 



a=l \fe=l 

sin(Ai + £ - z a ) sin(;e a + ^ + A i2 ) 
sin(Ai + £ + x a ) sin(x Q + 

Noting the relations (15.11) and using the same method as in [6], we find that if the 
parameters {x k } satisfy the following Bethe Ansatz equations 

2 ^ 2M j 



7 — 1 k — 1 



M 



2^ . . - s . , r, a = l,---,M, (5.7) 

f-^ sin i„ + x fc ) sm(x Q - x fc ) 
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the two sets of states • • • ,x M ) constitute the entire eigenstates of the Gaudin oper- 

ators 



HjV®(x ir --,x M ) = E j ¥ i \x l ,---,x M ) 



i = 



1,2. 



(5.8) 



The functions i£j is 



2 M 



sin 2zj 



Ej = cot 2zj + ^2 cot (^j + £ ~ Zj) + 



sin(x fc — Zj) sin(xfc + Zj) 



(5.9) 



j=l k=l 



6 Conclusions 

We have studied the second reference state of the open XXZ spin chain with non-diagonal 
boundary term, which leads to the second set of Bethe state (14. 8p . These Bethe states give 
rise to the corresponding Bethe Ansatz equations (14.101) and eigenvalues (14.111) proposed in 
[10. 9j by the functional Bethe Ansatz method. In the quasi-classical limit, two sets of Bethe 
states (15.51) and (I5.6P constitute the complete eigenstates of the associated Gaudin model. 

Very recently, an exact solution of the eigenvalue of the transfer matrix of open XXZ spin 
chain for arbitrary boundary parameters was proposed by functional Bethe Ansatz [18J and 
by representation of q-Onsager algebra [19] . It would be interesting to rederive their results 
in the framework of algebraic Bethe Ansatz. Moreover, such structure of multiply reference 
states found here also appears in open spin chains associated with higher rank algebras [20J. 
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